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ANNULAR  ELECTRON  BEAM 


Han  S.  Uhm 

Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  Md.  20910 
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A  self-consistent  theory  of  the  free  electron  laser  instability  Is 

developed  for  a  hollow  electron  beam  propagating  through  an  undulator 

(multiple  mirror)  magnetic  field.  The  stability  analysis  is  carried 

out  within  the  framework  of  the  linearized  Vlasov-Maxwell  equations. 

It  is  assumed  that  the  beam  Is  thin,  with  radial  thickness  much  smaller 

than  the  mean  beam  radius,  and  that  v/yj^  «  1,  where  v  Is  Budker’s 

2 

parameter  and  Yjj^c  is  the  characteristic  energy  of  the  electron  beam. 

The  dispersion  relation  describing  the  free  electron  laser  instability 
in  a  hollow  relativistic  electron  beam  is  obtained  for  an  equilibrium 
distribution  function  In  which  all  electrons  have  same  value  of  trans¬ 
verse  energy  and  the  same  value  of  canonical  angular  momentum,  and  a 
Lorentlan  distribution  In  axial  momentum.  It  is  shown  that  the  influence 
of  finite  radial  geometry  plays  a  critical  role  in  determining  detailed 
stability  behavior.  Moreover,  the  growth  rate  and  bandwidth  of  the 
instability  can  be  expressed  in  terms  of  Budker's  parameter  v.  Instead 
of  the  plasma  frequency  as  In  the  case  of  a  uniform  density  beam. 
Furthermore,  it  Is  found  that  free  electron  laser  stability  properties 
exhibit  a  sensitive  dependence  on  axial  momentum  spread. 


I.  INTRODUCTION 


In  recent  years,  there  has  been  a  growing  Interest  in  the  free 

electron  laser  instability^  ^  in  connection  with  Intense  radiation 
8*10 

generation.  For  the  most  part,  previous  theoretical  analyses  of 

this  instability  have  been  carried  out  for  an  electron  beam  with 
uniform  density,  neglecting  the  Influence  of  finite  radial  geometry. 
Although  this  is  a  reasonable  first  approach  to  the  problem,  for 
detailed  application  to  present  experiments  it  is  necessary  to  inves¬ 
tigate  the  important  Influence  of  finite  radial  geometry  on  stability 
properties.  In  this  regard,  in  the  present  analysis  we  investigate 
free  electron  laser  stability  properties  for  an  annular  electron 
beam  propagating  through  an  undulator  (multiple  mirror)  wiggler  field, 
including  the  full  Influence  of  finite  radial  geometry. 

Equilibrium  and  stability  properties  are  calculated  for  the  choice 
of  equilibrium  beam  distribution  function  [Eq.  (19)] 

fJ(H,  Pg,  p^)  -  K6(P0  -  Po)F(pf)G(p^), 

where  H  is  the  energy,  P^  is  the  canonical  angular  momentum,  p  = 

2  2  1/2 

(P  Pa)  transverse  momentum,  p  is  the  axial  momentum, 

r  0  .  * 

and  K  is  a  normalization  constant.  The  present  analysis  is  carried 
out  within  the  framework  of  the  linearized  Vlasov-Maxwell  equations, 
assuming  that  the  beam  thickness  is  much  less  than  the  mean  beam 
radius  Rq,  and  that  v/Yjj  1»  where  v  is  Budker's  parameter. 

It  is  also  assumed  that  the  electron  beam  propagates  through  an 

undulator  (multiple  mirror)  magnetic  field  and  that  the 

amplitude  B  of  the  wiggler  field  is  small  in  comparison  with  the  average 


field  Bq. 


The  basic  assumptions  and  equilibrium  properties  are  discussed 
In  Sec.  II.  The  formal  stability  analysis  for  azlmuthally  symmetric 


3 


perturbations  0/36  =  0)  Is  carried  out  In  Sec.  Ill  for  general 

2 

transverse  and  longitudinal  distribution  functions,  F(pj^)  and  G(p^). 
The  matrix  eigenvalue  equation  (65),  when  combined  with  Eq.  (67), 
constitutes  one  of  the  main  results  of  this  paper  and  can  be  used  to 
Investigate  stability  properties  for  a  broad  range  of  distribution 
functions  corresponding  to  a  thin  annular  electron  beam. 

In  Sec.  IV,  the  dispersion  relation  for  the  free  electron  laser 
Instability  Is  obtained  for  the  specific  choice  of  equilibrium 
distribution  function  in  which  all  electrons  have  the  same  value  of 
transverse  momentum,  and  a  Lorentzlan  distribution  in  axial  momentum 
(Eq.  (69)] 


G(P,) 


•  (  P,  -  P^)^  + 


Here,  A  Is  the  characteristic  axial  momentum  spread  about  the  mean 
value  p.  «Y.mV.  .  One  of  the  most  Important  consequences  of  finite 

D  D  D 

radial  geometry  Is  the  fact  that  the  maximum  Instability  growth  rate 
occurs  for  a  value  of  beam  radius  Rg  satisfying  [Eq.  (97)] 


®0t^“0s 


where  Is  the  radius  of  outer  conducting  wall,  Is  the  1th  root 

of  “Os  '^l^“0s^  “  Moreover, 

the  growth  rate  and  instability  bandwidth  can  be  expressed  in  terms 
of  Budker's  parameter  v  and  the  amplitude  of  the  wlggler  field  B.  It 
Is  found  that  the  maximum  growth  rate  Increases  with  increasing 
beam  Intensity  (v) . 


4 


The  dependence  of  stability  properties  on  the  axial  momentum 
spread  (A)  Is  also  investigated  In  Sec.  IV.  It  is  shown  that  the 
maximum  growth  rate  of  the  Instability  decreases  as  the  axial  momen¬ 
tum  spread  A  is  Increased.  However,  the  instability  bandwidth  increases 
with  increasing  A.  He  conclude  that  free  electron  laser  stability 
properties  exhibit  a  sensitive  dependence  on  axial  momentum  spread  A. 
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II.  EQUILIBRIUM  CONFIGURATION  AND  BASIC  ASSUMPTIONS 

A.  Basic  Assumptions 

The  present  analysis  assumes  an  Intense  annular  electron  beam 
with  characteristic  thickness  2a  and  mean  radius  Rq  propagating  In 
the  z-dlrectlon  (Fig.  1).  We  Introduce  a  cylindrical  polar  coor¬ 
dinate  system  (r,6,z)  with  z-axls  along  the  axis  of  synraetry;  r  Is 
the  radial  distance  from  the  axis  of  symmetry,  and  6  Is  the  polar 
angle  In  a  plane  perpendicular  to  the  z-axls.  To  make  the  analysis 
tractable,  the  following  simplifying  assumptions  are  made. 

(a)  The  thickness  of  the  annular  electron  beam  Is  much  smaller 
than  Its  mean  radius,  l.e., 

a/Rg  «  1.  (1) 

(b)  It  Is  assumed  that  the  beam  density  and  current  are  suffi¬ 
ciently  small  that  equilibrium  space  charge  effects  are  negligibly 
small,  and  the  equilibrium  self  magnetic  field  can  be  neglected  in 
comparison  with  the  applied  magnetic  field  B^.  That  Is,  we  approxl- 
mate: 

E°  -  0  -  B°  ,  (2) 

'\,a  *  ^  ' 

where  and  are  the  equilibrium  electric  and  magnetic  self  fields, 
respectively. 

(c)  Consistent  with  Assumption  (b).  It  is  assumed  that 

v/Y^  «  1, 

2  2  2 
where  Yj^nc  Is  the  characteristic  electron  energy,  and  v  *  M^^e  /me 


(3) 
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is  Budker's  parameter.  Here,  -e  is  the  electron  charge,  c  is  the 

speed  of  light  in  vacuo,  m  is  the  electron  rest  mass,  and  N.  = 

b 

2t\  /q  dr  r  n^  (r,2)  is  the  number  of  electrons  per  unit  axial  length. 

For  a  thin  annular  electron  beam  with  radial  thickness  2a  and  constant 
density  n^,  the  inequality  in  Eq.  (3)  can  be  expressed  in  the  equivalent 
form 


WpRQ/c^  «  Ro/a.  (A) 

2  2 

where  (21^  =  ATrn^e  is  the  relativistic  plasma  frequency-squared. 

(d)  It  is  assumed  that  the  electron  beam  propagates  through  an 
undulator  (multiple  mirror)  magnetic  field.  The  azimuthal  component 
of  vector  potential  A0(r,z)  for  the  equilibrium  magnetic  field  can 
be  expressed  as 


A^Cr.z)  »  1  Bq  ^r  -  f-  I^(k0r)cosko2^  .  (5) 

where  Xg  =  2Tr/kg  is  the  periodicity  length,  Ij^(x)  is  the  modified 

Bessel  function  of  the  first  kind  of  order  unity,  and  R  *  ^qz^*^  * 

kgZ  *  ±ir)/Bg^(r  =  0,  k^z  ■  0)  is  the  mirror  ratio.  Here,  the  equilibrium 

magnetic  field  components  are  given  by  B-  *  -9A?/9z  and  B.  *  r  ^(9/9r) (rA^) . 

Ur  b  Oz  V 

(e)  We  further  assume  that 


«  1 


(6) 


where  a  is  the  half-thickness  of  the  annular  beam.  Defining  the  charac¬ 
teristic  wlggler  amplitude  by 


B 


2 


R-1 

R+1 


> 


the  equilibrium  vector  potential  [Eq.  (5)]  can  be  expressed  in  the 


equivalent  form 


AjCr.z)  -  f 


i w  1  ^1^^ 

2V  ^"BQr  I^(k, 


o%>  ""’"o' 


From  Eq.  (7),  the  components  of  the  equilibrium  magnetic  field 
t  “  ®0rlr  ^Oziz 


Bor(r.z)  -  -  |^  A°(r.z) 


where 


1  - 

“  ■  2  *^0^0®  Ij^CICqRo) 


=:  -  B^sink^z, 


where 


®0z^’'’*^  “  r  I?  I*A°(r,z)) 


1  , 

'o  -  2  Vo  iTV?  “ 


coskgZ 


g  .  1  k  R  g  .  b  ^  I  (k  r  ) 

®z  2  Vo®  I^CkpEp)  ®0  R+1  ^O'Vo^’ 

Note  from  Eq.  (9)  that  Bq  is  equal  to  the  average  axial  magnetic  field 
(averaged  over  the  periodicity  length  lir/k^).  Note  also  that  the 

approximate  expressions  for  Bq^  and  Bq^  given  in  Eqs.  (8)  and  (9)  are 
valid  over  the  radial  extent  of  the  beam  (Rq  'a<r<RQ'fa)  within 
the  context  of  the  assumption  that  kQa  «  1  [Eq.  (6)] 
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(£}  From  Eq.  (9),  there  is  an  oscillatory  (wiggler)  contribution 

to  the  axial  magnetic  field  proportional  to  -  (B  /B^.)  cosk.z.  Throughout 

z  u  u 

the  present  analysis.  It  is  assumed  that  the  amplitude  of  the  wiggler 
field  Is  small,  l.e., 

A 

B 

«  1  .  (10) 

The  Inequality  in  Eq.  (10)  assures  that  the  axial  modulation  of  the 

equilibrium  beam  envelope  is  small  .  Equally  important, 

Eqs.  (1)  and  (10)  assure  that  the  axial  momentum  p  =  ymv  ,  where 

z  z 

2  22  1/2 

Y  =  (1  +  )^  /m  c  )  ,  is  a  good  approximate  single-particle  invariant 

to  lowest  order.  In  Appendix  A,  we  present  a  detailed  investigation 
of  the  electron  trajectories  in  the  equilibrium  fields  described 
by  Eqs.  (2),  (8),  and  (9),  assuming  that  the  axial  motion  is  non¬ 
resonant  with 


1,2  2  j 

Vs  * 


(11) 


Here  v^  =  p^/vm  is  the  axial  velocity  of  a  typical  beam  electron,  and 
“c  “  relativistic  electron  cyclotron  frequency. 

For  example,  in  the  strong-magnetic-field  limit  (w^  »  Kqv^),  it  is 
shown  in  Appendix  A  that  pz  =  const,  whenever  the  inequalities 
[Eqs.  (A. 20)  and  (A. 21)) 


2  Rq  Bq  v/)  ^  •  <12) 

and 

i  ^  j  1 »  (13) 

are  satisfied.  From  Eqs.  (1)  and  (10),  the  inequalities  in  Eqs.  (12) 
and  (13)  are  straightforward  to  satisfy  in  the  parameter  regimes 
of  experimental  interest. 


.  .. 
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(g)  Finallyi  in  the  stability  analysis  (Sec.  Ill),  it  is  assumed 
that  the  waves  are  close  to  resonance  with 

lo)  -  (k  +  nkQ)V^^|  «  kjjV^,  (0^,  (14) 

where  u  is  the  complex  oscillation  frequency,  k  +  nk^  is  the  excited 

3  0 

wavenumber,  and  =  /d  pv^f^  is  the  mean  axial  velocity  of  the 
electron  beam. 

In  summary,  the  present  equilibrium  and  stability  analysis  is 
carried  out  within  the  context  of  the  inequalities  in  Eqs.  (1),  (2), 

(4) ,  (6) ,  and  (10)  -  (14) .  Moreover,  the  vector  potential  for  the 
equilibrium  magnetic  field  is  approximated  by  Eq.  (7). 


B.  General  Equilibrium  Properties 


Within  the  context  of  Assumptions  (b) ,  (e),  and  (f)  in  Sec.  II. A, 
the  equilibrium  (3 /3 1  =  0)  beam  distribution  f^(r,2,]^)  is  generally  a 
function  of  the  single-particle  invariants'^  corresponding  to  particle 
energy , 


H 


/  2  4  ^ 
(m  c  + 


2  2.1/2 
c  Z  ) 


(15) 


canonical  angular  momentum, 

Pq  -  r[pg  -  I  Aj(r,*)l  ,  (16) 

and  axial  momentum 

Pz  (17) 

where  Ag(r,z)  is  defined  in  Eq.  (7),  and  -  (PytPg.p^)  is  the 
mechanical  momentum.  In  Eqs.  (15)  -  (17),  we  have  neglected 
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equilibrium  self-field  effects  and  assumed  chat  the  wiggler 
amplitude  B  is  sufficiently  small  [Eqs.  (10),  (12),  and  (13)]  that 
is  an  approximate  invariant.  In  this  regard,  we  also  conclude 
that  the  transverse  momentum 

/„2,  2  2  2  2.1/2 

Pj^  =  (H  /c  -  m  c  -  p  )  ,  (18) 

z 


is  also  an  approximate  invariant.  Here,  the  particle  velocity  v  is 

'X/ 

2  222  1/2 

related  to  the  momentum  ^  by  v  *  ]^/y®  where  y  =  H/mc  =  (1+^/mc  ) 
is  the  relativistic  mass  factor. 

In  Che  present  analysis,  we  consider  the  general  class  of  annular 
electron  beam  equilibria  described  by^^ 


f°(r,z,^)  -  K5(Pg  -  Pq)F(pJ)G(p^) 


(19) 


where  K  is  a  positive  normalization  constant,  P^.  »  const,  is  defined  by 


p  _ .  1  2 

^0  2  c  *^0  » 


(20) 


and  Che  transverse  and  longitudinal  momentum  distribution  functions 
are  normalized  according  to 


•00  ^00 

dpj  F(p|)  =  1  ,  dp  G(p  )  =  1. 
J  Q  J  —00 


(21) 


Defining  r^^  “  cpj_/eBQ,  and  carrying  out  some  straightforward 
algebraic  manipulation,^^  the  electron  distribution  function  in  Eq.  (19) 
can  be  expressed  in  the  equivalent  form 
N, 


"  2  2  2,,  2  2, ,1/2 
TT  [(r  -rj^)(r2-r  )] 


(22) 


[6($  -  $q)  +  6($  -  IT  +  Iq)]  , 


X 
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where  $  -  8in“^(pQ/Pj^) ,  and  (for  B^/Bp  «  1)  the  quantities  r^^  r^  and 
can  be  approximated  by 


rj^(*)  - 


II  -  |- cost,,.)  (l- 


2r, 


cosk^z 

0 


B«  — 0‘)  J 


1/2 


(l  -  |-  co.k„t 


rjU)  -  rjU)  + 


and 


I'  ''  /ft  \  ’ 

(l  -  |-  cosk^z) 

r2(l  -  |-  coskQz)  -  ^ 


(ji.  -  sin 


2rr, 


In  Eq.  (22),  ^(x)  is  the  Heaviside  step  function  defined  by 

1  ,  X  >  0  , 

0  ,  X  <  0  , 


«(x) 


)■ 

(23) 


(24) 


(25) 


(26) 


and  is  the  number  of  electrons  per  unit  axial  length  defined  by 

Nb  -  2it  I  drr  |  d|  J  dp^p^  j  dp^  f®(r,z,]^)  .  (27) 


Here,  is  the  radius  of  the  grounded  cylindrical  conducting  wall. 
Several  equilibrium  examples  of  the  general  foirm  given  by  Eq.  (22) 
have  been  discussed  previously  in  the  llterature^^  in  the  absence  of 
wiggler  field  (b*0). 
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III.  LINEARIZED  VLASOV-MAXWELL  EQUATIONS 
A.  General  Formulation 


In  this  section  I  we  develop  the  general  formalism  to  Investigate 
stability  properties  for  perturbations  about  the  class  of  beam  equilibria 
described  by  Eqs.  (19)  and  (22).  Expressing  the  beam  distribution 
function  as  f^(;jc,]j,t)  =  f^  +  6fj^(jc,]^,t) ,  and  neglecting  equilibrium 
self-field  effects  [Sec.  II. A],  the  linearized  Vlasov  equation 
can  be  expressed  as 


V  X  \ 


-e(- 


5A 


X 


V 


)•  - 


(28) 


where  ^  =  ]g/ym  is  the  particle  velocity,  the  equilibrium  magnetic 
field  components  ■  ®0r€r  ^  ®0z^z  given  in  Eqs.  (8)  and  (9), 
and  use  has  been  made  of  the  Lorentz  gauge  to  express  the  perturbed 
magnetic  field  as  6jg(;|c,t)  *  X,  ^  perturbed  electric 

field  as  £jg(;|c,t)  ■  -  V4<^(;^,t)  -  (1/c)  (3/3  t)6/^(;|t,t) ,  where 

•  5^(;^,t)  +  (l/c)(3/3t)6^(;|c,t)  ■  0.  The  corresponding  linearized 
Maxwell  equations  that  determine  the  evolution  of  £A  and  6^  self- 
conslstently  are  given  by 

and 

6^  =  -Aite  J  d^p  6f|j  ,  (30) 

where  the  perturbed  current  -e/dp  v  Sf^^  and  perturbed  charge  density 
2 

-e/d  p6£j^  are  calculated  from  Eq.  (28).  In  the  present  analysis,  we 
assume  azlmuthally  symmetric  perturbations  with 
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“b(Wt)  ■  «  )  -  ilp.  64(;5',t') 

«'  “r  ■>  W-'')]  9  0,. 

*  I - 'b<«  •«  )  - 

where  d;jc'/dt’  ■  and  djg'/dt'  -  -e^'  ]^®(;^')/c,  and  the  particle 

trajectories  the  equilibrium  field  configuration  satisfy 

the  "initial"  conditions  ;|c'(t’  -  t)  »  jc  and  ]^'(t'  »  t)  -  (Appendix  A). 

We  now  simplify  the  right-hand  side  of  Eq.  (35).  Without  loss 
of  generality,  the  perturbation  amplitudes  are  expanded  according  to 

6^(r,2)  =  ^(ji^(r)expIi(k+nkQ)z]  ,  (36) 

n 

etc.  The  particle  orbits  (;|c',j^')  are  calculated  in  detail  in 

Appendix  A  for  the  case  of  small  wlggler  amplitude  (B  /B.  «  1) 

z  0 

assuming  that  the  axial  and  cyclotron  motions  are  far  removed  from 
cyclotron  resonance  [Eq.  (11)],  l.e., 

t.2  2  ,  2 

u  z  c 

where  is  the  axial  velocity  and  a»c  “  eBp/ymc  is  the 

relativistic  cyclotron  frequency.  To  lowest  order,  the 
axial  motion  is  free-streaming  [Eq.  (A. 8)] 


(t-  -  t)  ,  (3; 

and  the  radial  orbit  [Eq.  (A. 13)]  and  azimuthal  orbit  [Eq.  (A. 17)] 
contain  oscillatory  contributions  proportional  to  cosm^T,  sinu^T, 
sin(kQZ  +  I^qV^t)  and  cos(kQZ  +  where  t  *  t'  -  t.  For  present 

purposes,  in  the  t*  -  integration  on  the  right-hand  side  of  Eq.  (35), 
we  retain  terms  proportional  to 

[o)  -  (k  +  nk-)p  /ym]  ^  , 
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and  assume  that  the  value  of  [u  -  (1c  +  nkQ)p^/Ym]  is  well  removed 

from  resonance  with  the  cyclotron  motion  and  the  beam  motion  [Eq.  (14)],  i.e.. 


|u  -  (k  -»•  nkQ)p^/Ym|  «  ,  kgV^^  . 


Within  the  context  of  Eqs.  (1)  and  (14),  Eq.  (35)  can  be  approximated 


by 


2  222  1/2 

where  y  ■  H/mc  -  (1  +  /m  c  )  .  Moreover,  within  the  context  of 

Eq.  (14),  only  those  contributions  to  v^  and  Vg  proportional  to  sln(kQZ  +  ^qV^t) 


and  cos(k^z  +  k^v  t)  are  retained,  i.e.,  on  the  right-hand  side  of 
u  u  z 

Eq.  (38),  we  retain  contributions  to  v'  and  v*  of  the  form 

r  0 

2 


i*)_ 


^r  2  ^0  B.  2  . 

0  u,  -  krtV 


krtV  sln(knZ  +  k-v  x), 
z  z  u  z  u  U  z 


O’z 


(39) 


2  .22 
“c  '^O'^z 

- -  cos(kQZ  +  k^v^x)  ,  (40) 

u  -  knV_  c 
c  0  z 

where  v  *  p  /ym,  u  >  eB^/ymc,  and  use  has  been  made  of  Eqs.  (A. 13) 

2  2  C  U 

and  (A.  17).  Finally,  since  the  oscillatory  modulation  of  the 

radial  orbit  Is  small-amplitude  [Eq.  (A. 13)],  we  approximate 

r*  s  r  In  the  arguments  of  the  perturbation  amplitudes  on  the  right- 

hand  side  of  Eq.  (38),  I.e.,  6^(r',z')  e  £^(r,z'),  5A  (r',z')  s  6A  (r,s'), 

etc. 


r'  -  4. 


2  B„ 


Substituting  Eqs.  (36),  (37),  (39),  and  (40)  on  the  right-hand 
side  of  Eq.  (38),  and  making  use  of  the  various  approximations 
enumerated  above,  we  find  after  some  straightforward  algebra  that 
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can  be  expressed  as 


exp[i(k  +  nkQ)z]  J  /  '^*  \ 

■  '  J  I.  -  (k  -b  nl.„)»,  -  r  *..„j 


2  ,22 
c  0  z 


0  u  -  k_v  c 
c  0  z 


iRfl  B  “c  f 

2  :2  2  Vz 

II  /it  —  V  .  V  V  '  '  . 


0  Uj‘  -  k‘v 
c  0  z 


(Al) 


where  the  function  A^,(]^,a),k)  is  defined  by 

af?  3fS 

(j^.w.k)  *  2[Yinu>  -  (k  +  n'kQ)p^]  — |  +  (k  +  n’k^)  j—  .  (42) 

3Pi 

Evidently,  (j^.oi.k)  Is  an  even  function  of  v^  -  Pj./Yn*  Therefore, 

It  follows  from  Eq.  (41)  that 


-e  j  d^p  v^  «fb(i^c,llj)  -  0  . 


(43) 


That  Is,  the  perturbed  radial  current  Is  equal  to  zero.  The  radial 
component  of  Eq.  (29)  then  gives 


•  1  «  2  ^ 

•f-  ~  ^  (rA  )  -  (k  +  nk^)  A  +  A  -  0  , 
3r  r  ar  r,n  0  r,n  2  r,n  ’ 

c 


(44) 


for  the  perturbation  amplitude  A  (r).  Without  loss  of  generality, 

r  ,n 

for  the  TE  mode  polarization,  we  assume 

A^  ^(r)  -  0  ,  (45) 

r  pu 

In  the  subsequent  analysis.  The  Lorentz  gauge  condition  In  Eq.  (34) 
then  gives  the  simple  relation  between  A  (r)  and  4  (r) 

^  2  I A  ^ 

(k  +  »k(,)A,  „(r)  -  ^  ♦„(!)  .  (46) 

Making  use  of  Eqs.  (45)  and  (46),  £fb()(,]^)  can  then  be  expressed  as 


exp[l(k  nkQ)g]  j  r  ^  ^ 

"  *  n  ***  ~  (*^  +  ^  (k  +  nkQ)cj  ^n 


^0  B  “c  ‘^0^«  ,,  .  .  ,1 

®0  J  -  kV  “c  ^  8»“-l  Vl^,n+1^/  • 

c  0  z 


B.  Elsznvzlue  Eauatlons 


The  linearized  Maxwell  equations  are  obtained  in  this  section 
within  the  context  of  Eqs.  (1)  and  (14).  For  convenience  In  the 
subsequent  analysis i  we  define  the  dimensionless  potentials  as 


[♦»('>  -  r  • 


A__(i)  -  -Sj  A  (t)  . 
me 

Moreover,  we  also  Introduce  the  dimensionless  parameters 


2  ,,2^2  * 
“c  ‘  Vb 


and  the  effective  susceptibility 

xi^\(u,k)  -  Ave^  I  rdr  [dp -  ,  (52) 

Jo  J  (w  -  (k  +  nkQ)v^] 

where  X^,(]^,u,k)  la  defined  in  Eq.  (42),  and  h  eB^/yi^mc  In  Eq.  (50), 

2 

In  Eq.  (51),  Is  the  relativistic  mass  factor  defined  by  - 

(1  -  V?/c^)~^.  Noting  that  the  azimuthal  velocity  at  t'  ■  t  is  given 

D 

by  [Eq»  (A. 17)) 

Vg  -  -(2cyjj/y)Aco8kQZ  ,  (53) 
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and  making  use  of  Eqs.  (42)  -  (52),  the  Maxwell  equations  (29)  and  (30) 
can  be  expressed  In  the  approximate  form, 


{7  It  It  +  7  - i'^n 


^b^O 


(54) 


6(r  -  Rq) 


{r  3r  '  ar  "  ^2  ^2  ■  “*^0  ■*■  }  ^n+l^*^^  Rq 

•  {0«  ■  ■"  •  <”> 

{r  ir  '  at  ■  “2  “2  ■  (I*  +  "■‘o  ■  ''o'  l\-l' 

'  r  c  i 


Hr  -  R  ) 
,(r) - = - 


pn.n^n 


L(2) 


.(2) 


^n,n-l^n-l  ^  ^n,n+l^n+l 


]}• 


(56) 


for  a  thin  annular  beam  satisfying  Eq.  (1). 

Since  the  right-hand  side  of  Eq.  (54)  vanishes  except  at  r  >  R^, 
Eq.  (54)  can  be  expressed  as 


(rlT’^lT'"Pn)^n^’^>  “  0  ’ 


(57) 


for  r  Rq.  In  Eq.  (57),  the  parameter  p^  is  defined  by 


2  2,2  ,,  x2 

p^  -  w  /c  -  (k  +  nkp) 


The  solution  to  Eq.  (54)  is  given  by 

®NQ(p„r)  ,  Rq  <  r  <  Rj,  , 


(58) 


♦nCr) 


[cJo(Pnr)  , 


0  <  r  <  Rn  . 


(59) 


where  J.(x)  and  N  (x)  are  Bessel  functions  of  the  first  and  second 

X»  Iw 
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kind,  respectively,  of  order  i,  end  Is  the  radius  of  the  conducting 
wall  (Fig.  1).  In  Eq.  (59),  the  constants  A,  B,  and  C  are  determined 
by  the  boundary  conditions  that  the  effective  potential  vanishes 

at  r  "  R^,  and  that  ^^(r)  Is  continuous  at  r  ■  Rq.  Multiplying 
Eq.  (54)  by  r  and  Integrating  from  RQ(l-e)  to  Rpd+e)  with  e 
we  obtain  after  some  straightforward  algebra, 

(60) 

where  the  dielectric  function  T  ,  ((>>,k)  is  defined  by 

n  gii 


Jn*  (Pn®c^^'^n'  ^^n^O^ 

and  use  has  been  made  of  the  Wronsklan  Identity, 


(61) 


(64) 


J^,  (x)dN^,  (x)/dx  -  (x)dJ^,  (x)/dx  -  2/irx  .  (62) 

Similarly,  making  use  of  Eqs.  (55)  and  (56),  we  obtain 

^l,n+l^“’*^^^n+l%^  “T*^^n,n*n^V  “  ^“^^,n-l\-l^*0^  ’^n,n+l^n+l^^)^ ^ ^ 

(63) 

id 

/  ^  a  V 

i  Eq.  (63),  Is  a  function  of  •  u^/c^  -  (k  +  nkQ  +  kQ)^ 

iflned  In  Eq.  (58).  From  Eqs.  (60),  (63),  and  (64),  we  obtain  the 

^  A  ^ 

itrlx  equation  relating  ^^^.^(Rq),  and  ♦jj(Rq) 

r  4.11,2  (2)  JL,2  (2)  JL,  (1)  \  A 

ri,a+l'*’2A  »Xn,„+i  »  2  A  \  M 

Vl'V  I-  »• 

'  ^  ^  '"0,1.  J\^  ♦n<*oy 


n+l^*0^ 
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which  is  similar  in  general  form  to  the  matrix  equation  obtalnec  by 
Davidson  et  al.^  for  the  free  electron  laser  Instability  in  a  uniform 
density  electron  beam.  The  dispersion  relation  can  be  obtained 
from  the  condition  that  the  determinant  of  the  matrix  in  Eq.  (65) 
vanish. 

Evidently,  from  Eq.  (65),  an  explicit  evaluation  of  the  susceptibility 
defined  in  Eq.  (52)  is  required  for  a  detailed  stability  analysis. 

XI )  u 

Substituting  Eqs.  (22)  and  (42)  into  Eq.  (52)  gives 


“  2vmc^  dp^  dp^|2[Ym(i)  -  (k  +  n'kQ)p^] 


3F_ 
^  2 
3Pj. 


+  (k  +  n'ko)F  II 


;} 


(66) 


(Yb/Y)- 


[m  -  (k  +  nkQ)v^l 


2  2 

where  v  *  Nbe  /me  is  Budker's  parameter  defined  in  Eq.  (3),  and  use 
has  been  made  of  the  identity 

rR 


•((r2  -  r)(r  -  r^^)] 
dr  r - -  - 


[(r^  - 


TT 

2  * 


In  obtaining  Eq.  (66),  the  p^  derivatives  of  ^Kr^  -  r)(r  -  r^^)) 

X  t(r^  -  r^)(r2  “  r^)]  and  «(♦  -  <)q)  +  6((j>  -  n  +  have 
been  neglected,  since  the  corresponding  corrections  are  of  order 
|„  -  <k  +  nk(,)V^|/|«^  +  kjV^I  «  1,  or  smaller.  Equation  (66)  can 
be  further  simplified  by  integrating  by  parts  and  making  use  of 

F(0)  -  F(«)  -  G(«)  -  G(-»)  -  0  . 

Within  the  context  of  Eqs.  (3)  and  (14),  we  obtain 

1“^  -  O'  +  ■"'o)0‘  +  '"‘k 

F(p1)G(p^) 


[u  -  (k  +  nkQ)p^/YiB]' 


X 


(67) 
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after  some  straightforward  algebra.  Equations  (65)  and  (67)  describe 

stability  behavior  for  general  choice  of  equilibrium  distribution  function 
2 

^(Px)G(p^)*  and  can  be  used  to  Investigate  stability  properties  for  a 
broad  range  of  system  parameters. 
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IV.  STABILITY  PROPERTIES 


A.  General  Dispersion  Relation 


As  an  example.  In  this  section  we  Investigate  free  electron  laser 
stability  properties  for  an  electron  beam  with  perpendicular  distribu¬ 
tion  function, 

F(Pi)  *  5(Pi  -  Pxo)»  (68) 

where  pj^o  ~  const.  Is  related  to  the  radial  thickness  of  the  beam.^^ 
Moreover,  In  order  to  Investigate  the  Influence  of  axial  momentum 
spread  on  the  free  electron  laser  Instability,  we  assume  an  axial 
momentum  distribution  function  of  the  form 


^^^2^  “  n  .  .22’ 

(Pg  -  Pfe)  +  A 

where  A  Is  the  characteristic  momentum  spread  about  the  mean 


(69) 


momentum  p .  =  .  We  further  assume  that  the  momentum  spread  A  Is 

l)  D  D 

small  In  comparison  with  the  directed  momentum  p^^,  i.e., 


A  «  P^. 


For  small  values  of6p  =p  -  P.  consistent  with  Eq.  (70), 


approximate  the  term  p  ^/Yi®  in  Eq.  (67)  by 


_ 2 

ym 


-  + 

y.m  3 
^  Yb® 


u  f,  ^  2,  2  2.1/2 

where  Y^  "■  (i  +  Pb'™  ^  ' 

Substituting  Eqs.  (68)  -  (71)  into  Eq.  (67),  we  obtain 

(j)  .  (0)  2v  -  (k  -I-  nkQ)(k  +  n’kQ)c^  ] 

^^0,0’  Yb  [u  -  (k  +  nkQ)V^^  +  i|k  +  nkglA/v^m]^ 


(70) 


we  can 


(71) 


(72) 
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within  the  context  of  Eq.  (14).  Equation  (65)  can  be  sliiq>llfled 
considerably  by  making  use  of  Eq.  (72).  The  condition  for  a  non¬ 
trivial  solution  to  Eq.  (65)  Is  that  the  determinant  of  the 
matrix  vanish.  This  gives  the  general  dispersion  relation 


which  can  be  used  to  determine  the  complex  elgenfrequency  u  In  terms  of 
k  +  nkQ,  kQ,  V,  A,  and  A .  What  Is  most  remarkable  Is  that  Eq.  (73) 

Is  very  similar  In  general  form  to  the  result  obtained  by  Davidson  et  al.^ 
for  a  uniform  density  electron  beam. 

Consistent  with  Eqs.  (3)  and  (14),  the  elgenfrequency  u  can  be 
approximated  by  u  *  (k  +  nkQ)Vj^.  We  therefore  approximate  p^  and 

Pn+1 


Pn  *  . 

and 

Pn+1  "  ^Vl  “  ^  2ko(k  +  nkp)  +  k2]l/2  ^ 


(74) 

(75) 


Making  use  of  Eqs.  (74)  and  (75),  the  dielectric  functions  r, 

1  ,n+i 

and  Fq  ^  can  be  expressed  as 


(76) 


and 


O.n  ■  2  '•<‘>nV  =  2  -  'o<<^*0>?o'''n''c) 


oW'hW 


(77) 


for  k  +  nkQ  >  0.  In  Eqs.  (76)  and  (77),  I|^(x)  and  Kj^(x)  are  modified 
Bessel  functions  of  the  first  and  second  kind,  respectively,  of  order  I, 
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and  h(q  R^)  Is  the  effective  electrostatic  wave  admittance.  Note  from 
^n  0 

Eqs.  (76)  and  (77)  that  the  dielectric  functions  ^  and  1^^ 
are  generally  nonzero  mono tonic  Increasing  functions  of  and 

q^^^RQ,  respectively.  In  this  regard,  for  small  wlggler  amplitude 
(A.  <<  1)  ,  we  Investigate  free  electron  laser  stability  properties 


for  u  and  k  +  nk^  near  the  simultaneous  zeroes  of  the  transverse 
dispersion  relation  F.,  =  0,  and  the  longitudinal  dispersion 

relation  F*  +  (Tt/2Yv)x^^^  ®  0*  Th®  general  dispersion  relation  in 

u }  ti  D  n  y  Hi 

Eq.  (73)  can  then  be  approximated  by  the  simplified  form 


r  fr  +  v^o>l 

l,n-l[  0,n  2  \i,nj 


IT  .2  (0)  „ 

oA  X 

2  n.n-1  0,n 


Equation  (78)  Is  the  form  of  the  dispersion  relation  used  In  the 
remainder  of  Sec.  IV. 


B.  Linear  Elgenmode  Properties 


For  very  small  wlggler  amplitude  (A  0) ,  or  for  very  low  beam 
density  (v  -*■  0) ,  it  Is  evident  from  Eq.  (78)  that  the  linear  dispersion 
relation  for  transverse  electromagnetic  perturbations  Is  given  by 


and  the  linear  dispersion  relation  for  longitudinal  perturbations  Is 


given  by 


Here,  we  assume  k  +  nk^  >  0,  and  the  elgenfrequency  u  satisfies 
u  -  (k  +  nkQ)V^  for  a  low-density  beam.  Equation  (79)  can  be  expressed 


In  the  equivalent  form 
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^  -  (k  +  nkp  -  kg)^  -  -|2.  ,  (8: 

c 

where  Og^  is  the  s'th  root  of  “  0,  and  use  has  been  made  of 

the  definition  in  Eq.  (61).  Equation  (81)  is  the  familiar  vacuum 

waveguide  electromagnetic  dispersion  relation.  For  future  reference, 

Taylor  expanding  the  electromagnetic  dielectric  function  F,  , 

i,n— X 

in  Eq.  (61)  about  [u^/c^  -  (k  +  nkg  -  kg)^]R^  =  “Qg»  it  is  straight¬ 


forward  to  show  that  F^^  can  be  approximated  by 


„  /  .  .2  “Osl  r‘^2^“0s^  I* 

^l.n-l<“-*^^  -  ~  4  1^2  ■  “>^0  -  *^0^  “  j^2  J  Iji(Oo3Ro/Rc)J  ' 


Moreover,  making  use  of  Eqs.  (72)  and  (77),  it  can  also  be  shown  that 
the  effective  longitudinal  dielectric  function  is  given  by 

1  4.  -J- 

XT  2  r 

2y^  'O.n 

222  (®3) 

2y  l(k  +  nkg)‘^c'^  -  u  ) 

*  1  Q  2  * 

Y^h  [(*)  -  (k  +  nkg)V^  +  l(k  +  nkg|A/Y^m] 

The  expression  for  the  electrostatic  wave  admittance  h(q^RQ)  in  Eq.  (77) 

can  be  simplified  in  several  limiting  cases.  These  include: 

2  2 

(a)  Long  wavelength  perturbations  with  q  Rg^<  i*  in  this  case, 
h(q^Rg)  can  be  approximated  by 

*’^‘*nV  “  * 

(b)  Large  conducting  wall  radius  (l.e.,  R^/Rg-^  «). 


In  this  case 
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2  2 

(c)  Short  wavelength  perturbations  with  I"  this  case, 

hCq^Rg)  can  be  approximated  by 

It  Is  Instructive  to  examine  the  longitudinal  linear  dispersion 
relation  In  Eq.  (80)  for  perturbations  with  long  axial  wavelength 
and  A  -►  0.  Keep  in  mind  that  A-*-  0  corresponds  to  a  cold  axial  distribu¬ 
tion  function  G(p^)  =  ^(p^  ~  Making  use  of  Eq.  (83)  in  the 

cold-beam  limit,  we  find  that  Eq.  (80)  can  be  approximated  by 

[(*)  -  (k  +  nkQ)V^]^  =  ^  [(k  +  nkQ)^c^  -  u>^]«.n(R^/RQ)  ,  (84) 

^b 

12 

which  is  identical  to  the  dispersion  relation  obtained  by  Briggs 
for  space-charge  waves  in  a  thin  annular  electron  beam. 

To  further  orient  the  reader  and  illustrate  the  Interaction 
wavenumbers  and  frequencies  characteristic  of  the  free  electron 
laser  Instability,  we  consider  a  tenuous  electron  beam  with  v  0  and 
A  0,  and  determine  the  simultaneous  solutions  to  Eqs.  (79)  and  (80). 
Evidently,  making  use  of  Eq.  (83),  the  longitudinal  dispersion  relation 
in  Eq.  (80)  reduces  to 

uo  =  (k  +  nkQ)V^  ,  (85) 

for  V  0  and  A  0.  Solving  Eqs.  (81)  and  (85)  simultaneously,  we 
obtain  the  interaction  frequency 

“  -  <1  ‘  11  -  (1  -  •  <“) 

where  and  u  represent  the  upshlfted  and  downshifted  frequencies, 
respectively.  In  the  limit  of  short  wiggler  wavelength  satisfying 

% 
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2  2  2  +  + 

>>  Oq^,  the  upshlfted  frequency  u  and  wavenumber  (k  +  nk^)  can 

be  approximated  by  the  familiar  results 

»♦-(!+  S^)YjVb  •  <"> 

(k  +  nku)+  -  (1  +  et)Yjk„  , 

obtained  in  free  electron  laser  theories  neglecting  geometric  effects. 

For  sufficiently  long  wiggler  wavelength,  however,  it  is  evident  from 
Sq.  (86)  that  geometric  effects  play  a  significant  role  in  determining 
the  characteristic  frequency  and  wavenumber. 

Finally,  for  small  but  finite  wiggler  amplitude  and  beam  density, 
we  make  use  of  Eqs.  (72),  (82),  and  (83),  to  simplify  the  full  dispersion. 


relation  in  Eq.  (78).  This  gives 


2  9  “n 

^  -  (k  +  nkQ  -  kjj) - ^ 


u)  -  (k  +  nkQ)Vjj  +  i 


|k  +  nkQ|A^2 


(2ko(k  +  "*^0^  ■  ■  _2* 

c 


-  -4aA‘ 


“o«1 


where  use  has  been  made  of  Bq.  (81)  to  simplify  the  right-hand  side  of 
Eq.  (88).  Equation  (88)  is  the  form  of  the  nonlinear  dispersion 
relation  used  in  the  numerical  studies  in  Sec.  IV. C. 


C.  Numerical  Analysis  of  the  Dispersion  Relation 
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Defining  the  normalized  Doppler-shifted  elgenfrequency  by 


fi  -  [u  -  (k  +  nkQ)Vjj]/kQC  , 


the  dispersion  relation  in  Eq.  (85)  can  be  expressed  as 


■v^  «  +  ,  ^  21  3  «  -  n) 

0  \  J 


+  C  =  0  , 


where 


C  -  2(k  +  nkQ)/kQ  -  (k  +  nkQ)^/k2Yj  -  i  -  aJ^/k^R^  ,  (92) 

n  =  (k  +  nkQ)2A^/Y®m^kQC^  +  (2v/Yjh)t2(k  +  nkQ)/kQ  -  dog/kgR^  -  1]  . 


£  V  A  "^0  V/^ON 

■"^bV  *^0  ■  ■  kjR2/i^J^c-^  ''2<“0s> 


.  (94) 


A  careful  examination  of  Eqs.  (91)  and  (94)  shows  that  the  maximum 
coupling  between  the  electromagnetic  and  electrostatic  modes  occurs  at 
the  value  of  Rq/R^,  satisfying 

where  the  prime  (')  denotes  dJ^(x)/dx.  Making  use  of  the  identity 
J' (x)  ■  “  J, (x)/x,  Eq.  (95)  can  be  expressed  in  the  equivalent 


In  this  context,  we  find  that  the  maximum  growth  rate  occurs  for  the  value 
of  Rq/Rj.  given  by 


“  ®0f^“0s  ' 


where  Is  the  £th  root  of  ®  satisfying 
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Figure  2  shows  a  plot  of  the  electrostatic  wave  admittance  h 
versus  normalized  axial  wavenumber  (k  +  nkn)/k.  for  y.  •  5,  k.R  -  S, 

U  U  D  U  C 

and  Kq/Rq  ~  ^Ol^'^Ol  ~  note  from  Fig.  2  that  the  wave  admittance 

h  is  proportional  to  the  axial  wavenumber  (k  +  ^or  wavenumbers 

satisfying  (k  +  nkQ)/kQ  >  S.  For  short  axial  wavelengths  with  (R^  -  R^) 

X  (k  +  nkg)  »  it  is  readily  shown  from  Eq.  (77)  that  the  wave 
admittance  h  can  be  approximated  by 


V  V 


For  free  electron  laser  applications  characterized  by  axial  wavenumber 

2 

k  +  nkg  =  (1  +  6b)Y|jkg,  Eq.  (98)  constitutes  an  excellent  approximation. 

Shown  in  Fig.  3  for  e  *  -0.1, and  in  Fig.  4  for  e  ■  +  0.1,  are  plots 
of  (a)  normalized  growth  rate  0^  *  ImO  and  (b)  Doppler-shifted  real 
frequency  Of  »  ReO  versus  the  normalized  axial  wavenumber  (k  +  nkg) /kg, 
obtained  from  Eq.  (91)  for  A  -  0  and  parameters  otherwise  identical  to 
Fig.  2.  In  Figs.  3(b)  and  4(b) »  is  plotted  only  for  the  range  of 
wavenumber  (k  +  nkg) /kg  corresponding  to  instability  (fi^^  >  0).  Evidently, 
from  Figs.  3(a)  and  4(a),  the  maximum  growth  rate  and  the  range  of 
(k  +  >^bg)  corresponding  to  instability,  increase  with  increasing 
value  of  beam  intensity.  It  is  also  evident  from  Figs.  3(a)  and 
4(a)  that  the  instability  occurs  at  values  of  axial  wavenumber  satisfying 
(k  +  nkg) /kg  <  2y^  for  e  <  0,  and  (k  +  nkg) /kg  >  for  e  >  0. 

Moreover,  the  normalized  Doppler-shifted  eigenfrequency  satisfies 
|(2|  «  1,  which  is  consistent  with  Eq.  (14). 

The  dependence  of  stability  properties  on  the  axial  momentum 
spread  (A)  is  Illustrated  in  Fig.  5(a)  for  t  *  -0.1,  and  in  Fig.  5(b) 
for  c  *  +  0.1,  where  the  normalized  growth  rate  0^  is  plotted  versus 
(k  +  nkg) /kg  for  v/yjj  •  0.05,  several  values  of  A/yj^mc,  and  parameters 
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otherwise  Identical  to  Fig.  2.  Several  features  are  noteworthy  from 

Fig.  5.  First,  the  maximum  growth  rate  decreases  as  the  axial  momentum 

spread  Is  Increased.  Evidently,  the  growth  rate  Is  reduced  substantially 
by  introducing  a  very  small  amount  of  momentum  spread.  From  Eq.  (91), 
this  feature  represents  a  general  tendency  for  all  radial  mode  numbers  s. 
Second,  the  Instability  bandwidth,  l.e.,  the  range  of  (k  +  nkQ)/kQ 
corresponding  to  instability.  Increases  as  A/vj^mc  is  Increased.  Third, 
for  a  specified  value  of  e,  the  maximum  growth  rate  occurs  at  approximately 
the  same  value  of  axial  wavenumber,  regardless  of  the  momentum  spread 
(e.g.,  at  k  +  nk^  =  45.3  k^  for  e  =  -0.1,  and  at  k  +  nkg  =  53  kg  for 

e  ■  0.1).  Finally,  we  also  note  that  a  sign  change  of  the  parameter  £ 

shifts  the  unstable  region  in  (k  +  nkQ)  space.  Moreover,  for  a  given 
value  of  A,  the  maximum  growth  rate  for  e  =  -0.1  is  slightly  larger 
than  that  for  c  0.1.  In  general,  we  conclude  that  free  electron 
laser  stability  properties  exhibit  a  sensitive  dependence  on  the  axial 
momentum  spread  A. 
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V.  CONCLUSIONS 

In  this  paper,  we  have  formulated  a  self-consistent  theory 
of  the  free  electron  laser  Instability  for  an  annular  beam  propagating 
through  an  undulator  magnetic  field.  The  stability  analysis  has  been 
carried  out  within  the  framework  of  the  Vlasov-Maxwell  equations.  The 
equilibrium  properties  and  basic  assumptions  were  summarized  in  Sec.  II, 
and  the  formal  stability  analysis  for  azlmuthally  symmetric  perturbations 
(9/96  -  0}  was  carried  out  in  Sec.  Ill,  including  the  Important  Influence 
of  finite  radial  geometi'y.  In  Sec.  IV,  the  dispersion  relation 
was  obtained  for  perturbations  about  an  annular  beam  equilibrium 
in  which  all  electrons  have  the  same  value  of  transverse  momentum 
and  the  same  value  of  canonical  angular  momentum,  and  a  Lorentzlan 
distribution  in  axial  momentum.  One  of  the  most  important  features  of 
the  analysis  is  that  the  maximum  growth  rate  occurs  for  *  ^Oi^^Os' 

Moreover,  the  growth  rate  and  instability  bandwidth  can  be  expressed  in 
terms  of  Budker's  parameter  v.  Instead  of  the  beam  plasma  frequency 
as  in  the  case  of  a  uniform  density  beam.  In  Sec.  IV,  it  was 
shown  that  stability  properties  exhibit  a  sensitive  dependence  on  axial 
momentum  spread  L. 
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Fig.  1 
Fig.  2 

Fig.  3 

Fig.  4 

Fig.  5 


FIGURE  CAPTIONS 


Equilibrium  configuration  and  coordinate  system. 

Plot  of  electrostatic  wave  admittance  h  versus  (k  +  nkQ)/kQ 
for  -  5,  kjjR^  -  5,  and  R^/R^  -  ®01^“oi  “ 

Plots  of  (a)  normalized  growth  rate  (2^  and  (b)  Doppler-shifted 
real  frequency  versus  (k  +  nkQ)/kQ  [Eq.  (91)]  for  e  -  -0.1, 
d  "  0  and  parameters  otherwise  Identical  to  Fig.  2. 

Plots  of  (a)  normalized  growth  rate  (2^  and  (b)  Doppler-shifted 
real  frequency  versus  (k  +  nkjj)/kQ  [Eq.  (91)]  for  e  ■  0.1, 

&  >  0,  and  parameters  otherwise  identical  to  Fig.  2. 

Plots  of  normalized  growth  rate  0^  versus  (k  -t-  nkQ)/k()  [Eq.  (91)] 
for  v/Yjj  *  0.05,  (a)  c  ■  -0.1,  (b)  e  ■  0.1,  several  values  of 
A/Vj^mc,  and  parameters  otherwise  Identical  to  Fig.  2. 


APPENDIX  A 


ELECTRON  TRAJECTORIES  IN  AN  UNDULATOR  MAGNETIC  FIELD 


In  this  appendix,  we  determine  the  particle  trajectories 

(r' ,  0',  z' ,  p' ,  Pq,  p*)  that  pass  through  the  phase-space  point 
r  o  z 

(r,  6,  z,  p^,  Pq,  p^)  at  time  t'  =  t  for  an  electron  moving  in  the 
undulator  wiggler  field  described  by  [Eqs.  (8)  and  (9)], 


r 


■  2 


(A.l) 


=  B-  -  B  cosk*z, 
z  0  z  0 


(A.2) 


where  A_  =  27r/k-  is  the  wiggler  wavelength,  B  =  const,  is  the 
u  u  z 

amplitude,  and  Bq  =  const,  is  the  average  axial  magnetic  field. 

2  22422  1/2 

Making  use  of  '  const.,  and  ymc  =  (n  c  +c  p'  ) 

=  const.,  the  exact  equations  of  motion  can  be  expressed  as 


r 


I 


d9' 

dt' 


2Ymc 


iiV') 


coskpz'^ 


(A.  4) 


d  z’ 


eB, 


dt 


,2 


2Ymc 


('■  !S)  S; 


sink„z' 


(A.  5) 


where  v'  »  p'/Ym  =  dr'/dt’,  vl  “  Po/Ym  =  r'd0'/dt',  v'  »  p’/Yo  “ 
r  r  y  V  z  z 

dz'/dt',  and  use  has  been  made  of  the  expression  for  Ag(r',  z')  in  Eq.  (7). 


! 


We  now  examine  solutions  to  Eqs.  (A. 3)  -  (A. 5)  that  pass  through 
the  phase  space  point  (r,6,z,p^,Pg,p^)  at  time  t'  *  t,  assuming  that  the 
wiggler  amplitude  Is  small  with 


-  1  k  R  I-  «  1 

Bq  2  0  0  Ij^CRqRq)  Bq 


(A.  6) 


Moreover,  in  the  subsequent  analysis  we  denote  the  small  radial  excur¬ 
sions  about  r'  ~  Rg  by  the  radial  coordinate 


P’  -  r*  -  R^. 


(A.  7) 


To  lowest  order  in  B/Bg,  we  find  from  Eq.  (A. 5)  that  the  axial  motion 
is  free  streaming  with 


(A.  8) 


p’  *•  p  . 
*^2  Z 


(A.  9) 


Substituting  Eq.  (A. 4)  into  Eq.  (A.3)  and  making  use  of  Eq.  (A. 7), 
the  radial  equation  of  motion  can  be  expressed  as 


df2  vW  [vp*  ' 


0+P’>  ^  Bg^O  Ig(kgRg) 


’  rJ  B  *^o  ^iVo'^V^ 

- ^ - T  +  1  +  -g-  p  r  I - p  \ -  cosk-z'  (A.  10 

(Rg+p')2  ®0VP  hW  0 


B_  .  p  . 

B«  •‘O^O  I,  (k„R«)  ‘^®**^0* 


f 
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Expanding  the  right-hand  side  of  Eq.  (A.  10)  for  B  «  and  |p'|  <<  R^, 
and  making  use  of  Eq.  (A. 8),  we  find  to  lowest  order 


2 

P 


1  A. 

2  c 


“o' 


k  —  (t'  - 
0  Y®  ^ 


t) 


(A. 11) 


where  o)^  =  eB^/Ymc  is  the  relativistic  cyclotron  frequency.  The 

radial  motion  determined  from  Eq.  (A. 11)  exhibits  a  strong  resonance 

whenever  k_v  =  u  ,  where  v  =  p  /y®  is  the  axial  velocity.  In  the 
u  z  c  z  z 

present  analysis,  we  assume  that  the  values  of  kQ,  and  the 

characteristic  axial  velocity  are  such  that  resonance  does  not  occur. 


l.e. , 


,2  2 
knV 

0  z 


(A.  12) 


Integrating  Eq.  (A. 11)  subject  to  the  boundary  conditions  (p'.pp  = 
(PiPj)  at  time  t'  =  t,  the  radial  motion  is  described  by 


r 

p'  =  pcosoi)  T  +  -  slno)  T 

c  Y®®  c 


2 

+  2  ^0  I:  2  V  2  Vz">  (A-  ^3) 

0  “c  ■  ‘^O’^z 

/  ‘"O'^z 

-  I  cosk«zcos(jo  T  + -  sink^zsinu)  t 

\  0  c  0  c 


p^  «  p^cosu^T  -  Y®®j.Pslnu)^T 

^2 

i  It  2  ~ ,  2  2  Vz^^ 

0  u  -  k_v 
c  0  z 

+  (u)  cosk-zslnu)  t  -  kf,v  slnk.zcosu  t)], 

C  w  C  V  Z  V  c 
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where  v  ■  P./y®  and  t  =  t'  -  t.  Note  from  Eq.  (A. 13)  that  the  radial 
z  z 

excursions  about  r  -  remain  small  amplitude  as  long  as 


(1) 


B _ 

Bn  2  .22 

°  “c  ■ 


2a 


"  *0  ’ 

where  2a  Is  the  characteristic  radial  thickness  of  the  annular  electron 


(A.  15) 


beam.  Finally,  for  the  6 '-motion,  we  expand  Eq.  (A. 4)  for  B  «  B^ 
and  |p' I  «  Rq.  To  lowest  order,  this  gives 

de'  ^  1  B 

B  J.i  —  T-  /A  ■ 


2  “c  B^  • 


df  "c  Rq 

Substituting  Eq.  (A. 13)  Into  Eq.  (A.16),  we  obtain 

de*  p  ^r  . 

-T~r  ”  n  w  cosfa)  T  +  — sintij  t 
dt  Bq  ^  Y®Bq  c 


(A.16) 


u 


1  B _ 

2  B-  2  ,22 

°  “c  ■  ^0^2 


U  COSk^ZCOSU)  T 

c  0  c 


1,2  2 

+  k«v  sink-zsinoj  t - cosCk^z  +  k, 

0  z  0  c  u)  0 


o^’>] 


(A. 17) 


In  addition,  the  azimuthal  momentum  Is  readily  determined  from  p'  > 
ymr'de'/dt'  =  yvxa^p'  -  (Ym(D^RQ/2) (B/BQ)co8kQz’ ,  where  p'  and  z'  are 
defined  In  Eqs.  (A. 13)  and  (A. 8),  respectively. 

We  conclude  this  appendix  by  noting  that  Eq.  (A. 5)  can  be  used 
to  calculate  corrections  to  the  lowest-order  z-motion  In  Eqs.  (A. 8)  and 
(A.9).  Expanding  the  right-hand  side  of  Eq.  (A. 5)  for  |p'|  «  Rq 

A 

and  B  <<  B^,  and  making  use  of  Eq.  (A.16)  gives 


» 


(A. 18) 


where  p'  =  Y®dz'/dt’.  Substituting  the  expressions  for  z'  [Eq.  (A. 8)] 
z 

and  p’  [Eq.  (A. 13)]  Into  the  right-hand  side  of  Eq.  (A. 18)  and  Integrating 
with  respect  to  t',  we  find 

P;  '  Pz  a.^(kQRQ)(Ymu.^RQ)  |- 


1  B_ _ .  _  \ _ 1 _ 

R.  2  B-  2  .22  1  2  2  2 

(V  0  0  01^  -  /  0)^  -  IcqV^ 

X  |(i)^^cos  (RqZ  +  kQV^T)cosa)^T  -  cosk^z^ 


+  kQV^sin(kQZ  +  kQV^T)sin(D^  T 


f  Pr  1  B  *^cVz  ,  .  \  1 

I  Yinoj  Rf,  2  B„  2  .22  2  2 

|a)^^sln(kQZ  +  kQV^T)coS{u^T  -  slnkpZ^ 


(A.  19) 


1  B  ^O^z  I 

8  bT  ~2 - TT  tcos(2kQZ  +  2kQV^T)  -  cos2k  z]  )  , 

°  “c  -  ‘^O'^z  j 

where  t  ■  t*  -  t,  v  *  p  /ym,  and  p’  *  p  at  time  t'  =  t. 

Z  2  2  2 

Note  from  Eq.  (A. 19)  that  the  oscillatory  corrections  to  p'  =  p 

z  z 

[Eq.  (A. 9)1  are  small  whenever  the  wlggler  amplitude  B  is  sufficiently 

small  In  comparison  with  Bq.  Expressing  Eq.  (A.  19)  as  "*■  ^P^* 

the  Inequality  |fip^|  «  [p^l  is  hardest  to  satisfy  in  the  strong 

2  2  2  2  2  2 

magnetic  field  limit  where  u  >>  kAv  .  For  u  »  k»v  ,  we  find  from 

C  ^  Z  C  U  2 

Eq.  (A. 19)  that  t6p*|  <<  |p  |  provided  the  Inequalities 

Z  2 ' 
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and 


0)  R_ 
c  0 


) 


« 


1  , 


(A. 20) 


®o  * 

are  satisfied.  In  obtaining  £q.  (A. 20),  we  have  estimated  the 

characteristic  value  of  p  by  p  a,  where  2a  is  the  thickness  of 

the  annular  electron  bean.  The  inequalities  in  Eqs.  (A. 20)  and 

(A. 21)  are  straightforward  to  satisfy  for  experimental  parameters 

of  Interest.  Note  also  that  p*  =  p  whenever  Eqs.  (A. 20)  and  (A. 21) 

z  z 

are  satisfied,  and  the  axial  momentum  p  is  a  good  (albeit  approximate) 

z 

invariant . 


(k+nko)/k 


(k  +nko)/  k 


=  0.05 


